Localized charged states and phase separation near second order phase transition 
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Localized charged states and phase segregation are described in the framework of the phenomeno- 
logical Ginzburg-Landau theory of phase transitions. The Coulomb interactions determines the 
charge distribution and the characteristic length of the phase separated states. The phase separa- 
tion with charge segregation becomes possible because of the large dielectric constant and the small 
density of extra charge in the range of charge localization. The phase diagram is calculated and 
the energy gain of the phase separated state is estimated. The role of the Coulomb interaction is 
elucidated. 

PACS numbers; 



There is a common belief that complex interactions be- 
tween charge, spin, orbital, and lattice degrees of freedom 
may lead to inhomogeneous ground state with charge 
and phase separation. The role of these inhomogeneous 
states in the anomalous transport and magnetic prop- 
erties in manganites [H [3, S, 0, JE ^md in high tem- 
perature superconductors [7], |8|, |9|, [lOj is often discussed 
in literature. It was shown that inhomogeneous states 
may appear above Curie temperature and new charac- 
teristic temperature T* was introducedjl, Q. Note, that 
the tendency to the phase separation is widely discussed 
for manganitesjl, JR 11], as well as for high-temperature 
superconductors [7|, microscopic origin of this anoma- 
lous behavior is far from understanding. The origin and 
the temperature raiige of this phenomena is still an open 
problem [3, i, i, i, B, B 113, El. In this paper we dis- 
cuss the tendency and conditions of the formation of in- 
homogeneous states with the spacial charge localization 
and the phase separation within the phenomenological 
Ginzburg-Landau formalism and clarify the role of the 
Coulomb interaction in this phenomena. 

The problem of spatially inhomogeneous states in 
the charged systems frustrated by the Coulomb inter- 
actions is constantly under debate for different systems 
i, i, S, i, 0, El, m d, [H, 113, d. Usually the 
first order phase transition is considered where the charge 
density is coupled linearly to the order parameter (as an 
external field) [itI] or to the square of the order param- 
eter (local temperature) [l^ . Very often the Coulomb 
interaction is discussed on the qualitative level or on the 
last stage of the calculations [1, [H, [13, [l^ . The general 
method of consideration of the long-range Coulomb in- 
teraction for the two dimensional case was proposed in 
ref. [13, Ell . In the recent paper^ the effect of the long- 
range Coulomb interaction to the second order phase 
transition was considered. However the Coulomb interac- 
tion was considered only for the special cases, where the 
density distribution was approximated in the attempts to 
minimize Coulomb contribution to the free energy flO^, 
rather then it was obtained from the general formula- 
tion of the problem. In our opinion, The contribution of 
the Coulomb interaction is most important in the deter- 



mination of the form of the inhomogeneous states and 
therefore wc consider this interaction exactly. 

Phenomenological approach to the theory of phase 
transition, proposed by V.L. Ginzburg and L.D. Landau 
indicates that the properties of the system near the phase 
transition are determined by the closeness to the phase 
transition point and weakly dependent on the other prop- 
erties of the system. Therefore our results demonstrates 
that the properties of the spatially inhomogeneous states 
are determined by the closeness to the phase transition 
point, but not by the properties of the interactions in 
the system. We obtain the contribution of the Coulomb 
interaction to the thermodynamic potential of the inho- 
mogeneous state and show that this contribution may be 
relatively small, to allow the phase separation. Note, that 
the phase separation is possible for the system displaying 
phase transition of the second order. 

Let us consider doped system in the vicinity of the sec- 
ond order phase transition. We assume that the average 
concentration of the free carriers p is proportional to the 
dopant concentration x. The thermodynamic potential 
^ — J <j>{Vi P)dV describes the behavior of the order pa- 
rameter rj near second order phase transition and interac- 
tion of the order parameter with the charge density. Close 
to the phase transition it has the form: 
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where is the density of the thermodynamic potential 
in the high-temperature phase, 0,, is the density of the 
thermodynamic potential in the low-temperature phase, 
a, /?, D are coefficients in the expansion of the ther- 
modynamic potentials in powers of the order parame- 
ter [a — a (T — Tc), where Tc critical temperature in 
the absence of doping, a — 1/C, where C is the Curie 
constant), D is proportional to the diffusion coefficient. 
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Density of the thermodynamic potential (jjint describes 
the interaction of the order parameter with the charge, 
(Ti and (72 are the constants of interactions (the term 
C2V^P^ is important to provide global stability of the 
system[l3l, 1T2 > 0). The density of the thermodynamic 
potential (pcoui describes charging effects due to long- 
range Coulomb forces, p is the average concentration of 
the charge. We assume that Coulomb contribution is the 
strongest. The effect of average dopant concentration p 
to the thermodynamic potential in the high-temperature 
phase is included in 0o . The coefficient 7 is inversely pro- 
portional to static dielectric constant £, 7 = 1/e. As a 
result the coefficients in the thermodynamic potential de- 
pend on charge density p and the term aip{l + p(J2/(yi) 
determines the shift of the critical temperature of the 
phase transition due to local charge density. Therefore 
the influence of the electronic system on the phase tran- 
sition manifests itself as a continues shift of the critical 
temperature Tcp with the change of the local charge den- 
sity: 



Tcp(p) = Tc - (o-i/o(l + cF2p/(yi))/a 



(2) 



The local charge density may be different from p in some 
region of a finite size. 

In the future we perform analysis for 3D case. The 
analysis for 2D case is similar. Equilibrium state of the 
system is determined by two equations: d4>{ri, p)/drj — 0, 
d(j){ri, p)/dp = 0. The second equation determines the 
dependence of the charge density po on 77: 

Poiv) - P+^(l + 2a2p/ai)VV 

+ 2a2P/<Jimw? + v'^'v'm 



Substituting Eq.(3) to Eq.(l) we obtain the expression of 
the thermodynamic potential (1) (ppo for the equilibrium 
distribution of the charge density po{ri): 
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The negative sign in the third term of Eq.(4) (— (V?7^)^) 
indicates that uniform state may be unstable towards 
inhomogeneous fluctuations. Similar instabilities were 
discussed in the case of first order phase transition in 
Refs. [13,111. 

This instability leads to spatially inhomogeneous so- 
lutions. Most simple solution is spherically symmet- 
ric. This equilibrium inhomogeneous state, obtained by 
computer minimization of the thermodynamic potential 
Eq.(l), has spherical form and plotted in Fig.(l). This 
results are obtained by applying the numerical technique 
described in Ref.jl^]. This solution has characteristic 
size i?o and characteristic charge density inside of this 
region po under condition that average charge density in 
the system is equal to p. It is clear that the charge is 
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FIG. 1; (a) Spatial variations of the order parameter r\[r) 
and charge density A/9(r) = p(r) — p along the diagonal of 
the sphere in the equilibrium nonhomogeneous state, obtained 
by numerical minimization of the free energy; Qs is integral 
concentration of the charge as a function of distance from the 
center of sphere {Qs = Ap(r)S'r) 



concentrated near the surface of the sphere. It allows to 
use approximation of the double electrical layer for cal- 
culation of the Coulomb energy. At large distances from 
the sphere the order parameter and the charge density 
are equal to their equilibrium values p — p and rj — r]s (p) 
{ijs — or = —a{p,T)/P). As a result the density of 
the thermodynamic potential (t>s for this solution we can 
write in the form: 
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K{xo) 



Rt (5) 



Here we use dimensionless concentration x instead of 
charge density p {x = pa^ je, a is the lattice constant, 
and e is elementary charge) and x corresponds to the 
level of doping. The coefficients in this formula are de- 
fined as follows: A{xo) = ■n{a{xo,T))'^ / l3Vo, B{xo) ~ 
8TTD{-a{xo,T))/df3Vo, K{xq) ~ -fe'^dF{x - xq)'^ /a^Vo 
{a{xo,T) = a(T) + dix + dix'^, di = eai/a^, <72 = 
e^a2/aP, F is dimensionless factor which determines 
charge distribution F = 4a^ / <j)couidV / ^e^ R^dix — xq)^, 
F = 1/18 in the limit of the double electrical layer, Vq is 
the elementary layer that contains one sphere (Vq = V/n 
where n is the number of spheres), d is the thickness 
of the interphase boundary. Note, that Eq.(5) is valid 
if characteristic dimension of the nano-regions is larger 
then characteristic length ^ = hv/kBTcp{xo), where h is 
the Plank's constant, ks is the Boltsmann constant, v 
is the characteristic velocity of the problem (in the case 
of structural phase transitions it is equal to the sound 
velocity and in the case of electronic phase transitions it 
is equal the Fermi velocity). This may be achieved rela- 
tively easy if Tcp{xo) > lOK. Here we take into account 
the strong screening of the localized charges and write 
the energy in the double electrical layer approximation. 
In that case it is proportional to dRg. If the screening is 
absent this energy is proportional to i?g . The param- 
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etcrs A(xq) and B{xo) depend on temperature T, and 
parameter K{xo) depends on the average charge density 
X. Parameters {x, T) are external and we find the phase 
diagram in the space of these parameters. Since x is de- 
termined by the level of doping, we study the change of 
the properties of the system with doping. 

The conditions for the minimum of the potential 
4>s{Ro,xo) as a function of i?o and xq and the condition 
that minimum of the potential is lower then the potential 
of the uniform state (p^ < define the following set of 
equations: 



,iRo,xo)/dRo = 
s{Ro,xo)/dxo = 



= vo 



(6) 
(7) 
(8) 



These equations define the upper boundary of the ex- 
istence of the inhomogeneous state. Eq.(6) defines the 
equilibrium size of the charged domain Rq: 



Rn = 




2K{xo) 



(9) 



On the upper boundary of the region of the inhomoge- 
neous states (Eq.(8)) we obtain the localization radius 

Ro = 2A{xo)/3K{xn). As a result Eq.(8) gives the fol- 
lowing equation for xq as a function of temperature T: 

a{xo, T){ar + 2a2Xof + ^^^^'^^^ = o (10) 

Using Eqs.(lO) and (7) we write equation for upper 
boundary of the inhomogeneous states in the space of 
external parameters {x, T) (the sign ± is determined by 
the sign of cii): 



X = Xo{T) ± ( 



na 



36je^/3FD 



Y'^a^'^{xa{T),T) = (i (11) 



These equations define the temperature Tps\ of the 
transition to the region of stable inhomogeneous states 
in the space of external parameters {x, T) {Tp^i = fi{x)) 
where the energy of this inhomogeneous state is lower 
then the energy of homogeneous state. Similar consider- 
ation leads to the equation which defines lower boundary 
of the inhomogeneous states Tps2 = f2{x). Note that 
lower boundary Tps2 will be always close to the temper- 
ature Tcp{x), because we do not gain, but contrary loos 
in energy 0^ when the bubble of the high-temperature 
phase is formed {(jiniv = 0) = 0, and (jinirio) < 0). It 
leads to the considerable difference between the forma- 
tion of the bubble with the low-temperature phase in 
the high-temperature region and the bubble of the high- 
temperature phase in the low-temperature region. The 
first one is energetically favorable and therefore the re- 
gion of the existence of this bubbles is considerably larger. 
It is important to note that Eqs.(10,ll) which determine 
the phase diagram do not depend on d. This fact allows 
us to avoid optimization of our analysis with respect to d. 
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FIG. 2: Phase diagram of the system with Tc = 350K, 
ai/a = 50K, aija. = 510^K, -ye'^fSFD/Tra'^a = 2.33T^: I 
is is the region of the stable high-temperature phase; II is the 
region of the stable high-tcmperaturc phase with metastable 
bubbles of the low-temperature phase; III is the region of 
the stable inhomogeneous phase and metastable homogeneous 
high-temperature phase; IV is the region of the stable inho- 
mogeneous phase and metastable uniform low-temperature 
phase; V is the region of the stable uniform low-temperature 
phase. 



Typical phase diagram of the inhomogeneous state is pre- 
sented in Fig. 2 for the case of of positive tri > 0. Phase 
transition to the nonhomogencous state represents typi- 
cal phase transition of the first order. Metastable inho- 
mogeneous phase appears at the temperature Tp^o which 
is much higher then the temperature of phase transition 
and it is shown in the phase diagram by dashed line. 
This line is determined from the conditions Eqs.(6,7) and 
A{xo)'^ = 4:K{xo,x)B{xq,x) instead of Eq.(8). The char- 
acteristic size of the charged regions in that case is de- 
termined by the condition R.q = A(xq)/2K{xq, x). 

Very often the Coulomb energy is estimated as e^/a ~ 
4eV, which corresponds to two elementary charges local- 
ized in neighboring unite cells. However in our case to 
make the phase separation possible it is sufficient to have 
only small variation of the charge density per unite cell 
in comparison with average charge density a^{pa — p) 
0.1 — 0.2e. Therefore we can estimate the Coulomb con- 
tribution to the free energy Up as well as the energy gain 
due to formation of the low-temperature phase (the 
second and the forth terms in Eq.(5)) per one unite cell 
in the nonhomogencous phase Ui = 3|(/)i|a"^Vo/47r_RQ. We 
assume that d is small and neglect surface term in Eq.(5). 
For the Coulomb contribution we obtain: 



^ 3e'^{xo-xf \ dRoF 



(12) 



The characteristic energy in the brackets is of the or- 
der Mo = 0.0003 — O.OOley, here we use static dielec- 
tric constant £ = 30 — 40, which is consistent with the 
measurements on manganitcs and high-temperature su- 
perconductors. The increase of the bubble size leads to 
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the increase of the energy Up by the order of magni- 
tude or more in comparison with uq, but on the other 
hand geometrical factor F decreases Up. {F = 1/18 in 
the hmiting case of the double electrical layer). There- 
fore Up is evaluated as 0.003 — O.OleV. The corre- 
sponding estimate for the energy gain due to forma- 
tion of the low-temperature phase is ~ 3fcsTc/4 
K - a3(a(xo,r = 0))V4/3 ~ a^T,riixo,T = 0)/C, 
where C = N{r]o{T = 0)a?'f /^ks and N = a^^. There- 
fore Ur, ^ 0.02 - 0.05ey for ~ 270 - 700i^. Therefore 
in our case the phase separation becomes possible. The 
typical length scale of the nano-regions is estimated as 
i?o - 1 - 20nm. 

The analysis of the pair distribution function obtained 
by neutron scattering shows that the charge density in 
manganites is localized on the scale of 3 to 4 interatomic 
distances. This is the case even if holes are localized in 
the form of polaronsf^, [2l|. It confirms that holes are 
localized not in the form of well separated polarons but 
in the form of nanoregions with finite charge density of 
polarons. The extra charge in that case is not more then 
0.1 — 0.2e per unite cell. This is in consistence with our 
estimates. This state is characterized as the state with 
nano-dimensional charge and phase separation. Dynam- 
ics of these charged nano-regions may lead to high value 
of the dielectric constant in the low frequency range [ll[ . 

Therefore we have shown that the second order phase 
transition with critical temperature equal to zero at 
the certain doping (quantum critical point) is unsta- 
ble with respect to formation of the spatially inhomo- 
geneous charged states. Within the phenomenological 
Landau theory we have shown that these states appears 
at some temperature Tpsi, which is substantially higher 



then the temperature Tcp (Fig. 2). As a result the phase 
transition becomes effectively first order phase transition. 
Note that the Coulomb interaction determines the charge 
distribution, the screening and the characteristic length 
scale of the nonhomogeneous states. The spatially inho- 
mogeneous state becomes possible in the systems with 
the large dielectric constants (30-40) and with relatively 
small charge density variations. The nonhomogeneous 
states were discussed also in the case of the phase transi- 
tions of the first order [13, [3. However in that case spa- 
tial variations are related to the existence of two minima 
in the free energy in the region of hysteresis of the first 
order phase transition. In our case the nonhomogeneous 
states are related to the spatial charge redistribution. 

In conclusion we note that the localized charged states 
and the phase separation appears even in the case of 
the second order phase transition. Properties of these 
states are described within the phenomenological theory 
of the phase transitions. The Coulomb interaction de- 
termines the spatial charge distribution, the screening 
and the characteristic length of charge localization. The 
charge and the phase separation becomes possible be- 
cause of large dielectric constants and relatively small 
spatial variation of the charge density. We derive phase 
diagram and estimate different contributions to the free 
energy. 
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